arXiv:hep-ph/9502374vl 23 Feb 1995 


BARI-TH/95-191 
UCT-TP-222/95 
January 1995 


Tests of factorization in color suppressed 
nonleptonic B decay modes 


P. Colangelo®, C.A. Dominguez^f], N. Paver^Q 


“ Istituto Nazionale di Fisica Nucleare, Sezione di Bari, Italy 

^ Institute of Theoretical Physics and Astrophysics, University of Cape Town, Rondebosch 
7700, South Africa 

^ Dipartimento di Fisica Teorica, Universita di Trieste, Italy, and Istituto Nazionale di Fisica 
Nucleare, Sezione di Trieste, Italy 


Abstract 


We propose a number of tests of factorization for the color-suppressed nonleptonic 
decay channels: B —> K^*') rj^ and B —> 77 ^. The relevant leptonic constants 

and form factors are calculated within the QCD sum rules approach. 


^John Simon Guggenheim Fellow 1994-1995 

^Also supported by MURST, Ministero dell’Universita e Ricerca Scientifica e Tecnologica, Italy 



1. Introduction 


It has been recently observed [Q, ^ that the available experimental data on the lon¬ 
gitudinal polarization of the hnal particles in the decay B K* J/\b; pi = T{B ^ 
K* J /'^)ll/^{B K* J/tb) = 0.84 ± 0.06 ± 0.08, together with the ratio Rj/^a = 

r(S —i> K* J/^>)/T{B K J/4/) = 1.64 ± 0.34, can severely constrain the semilep- 
tonic B —> form factors if these nonleptonic transitions are evaluated within the 

factorization approach. 

Neglecting penguin operators, the relevant effective Hamiltonian density is 

Tiw = Vcs^ch (Cl(/i) Oi + C 2 (/i) O2) + h.c. , (1) 

where Oi = {siCi)v-A{cjbj)v-A and O 2 = {siCj)v-A{cjbi)v-A (hj are color indices, V — A 
refers to 7/^(1 — 75), Gp is the Fermi constant and Vim are Cabibbo-Kobayashi-Maskawa 
matrix elements). The values of the Wilson coefficients at the next-to-leading order in 
QCD can be taken as ci(mf,) = 1.08 ± 0.01 and C 2 {mb) = —0.19 ± 0.02 [|, 

The above transitions are color-suppressed. After Fierz reordering Eq. o. in the frame¬ 
work of factorization and vacuum saturation dominance, the factorized amplitude is 

A{B ^ aW J/T) = -^v:,Vba2 < K^*^\s-f^{l-^,)b\B > < J/T|c7^(l-75)c|0 >, (2) 

where 02 = C2 + Ci/Nc, with W is the number of colors. Eq. @ involves the matrix 
element 

< e)|c7/.c|0 >= /j/vp Mj/q> e* , (3) 

where the leptonic constant fj/q, is directly obtained from the measured decay J/T — 
Equation (^) also involves the semileptonic B —> RV matrix elements {q = p—p')'^ 

/u2 _ t\j2 

< K{p')\s-f,b\B{p) >= {p + p'),F^{q^) + ^-J[Fo(g2) - F^{q^)]q,, (4) 

and 

^We adopt the Bauer-Stech-Wirbel parameterization where: Fo(0) = A(0), ^ 3 ( 0 ) = To(0) and 
AAq^) = [{Mb + MKA/^MK^Aiiq^) - [{Mb - MKA/2MK^]A2{q^). 
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< K*{p\ri)\s-f^,{l - -f^)h\B{p) > = 


2 V{q^ 


Mb “ 1 “ Mx* 


{Mb + MK*)Ai{q^)p 


2^^* 


A2{q^ 


■{r]* ■p){p + p)^ 


2M. 


{V* ■ p)=^q,{As{<l") - Mq^)) 


Mb + Mx* ' q 

Actually, in Eq. the form factors must be evaluated at = Mj^^ and, in particular, 
the quantities pi and-Rj/^ depend on the ratios: E(Mj/^)/Ai(Mjy,j), A2(Mj^,j,)/Ai(Mj/^) 
and F,{Ml^)/A,{Ml^). 


( 5 ) 


Using predictions of popular theoretical models for the form factors,^ difficulties are met 
in simultaneously reproducing the experimental outcomes for pi and Rj/q, B0. This 
has prompted a number of interesting phenomenological analyses trying to elucidate the 
source of the problem. On the one side, the blame could be put on factorization which, 
in principle, can receive substantial corrections from the nonfactorizable component of 
Eq. that has the large coefficient Ci 0 . At present, such effects are quite difficult to 
predict by a dynamical calculation. 


On the other side, the alternative attitude would be to accept factorization, which is an 
appealing exact feature of multicolor chromodynamics in the limit Nc ^ oo 0.0 and 
focus, instead, on form factor predictions. Actually, not all form factor models are quite 
predictive in the region of large ~ ^J^cl in these cases the g^-dependence has to 

be assumed, e.g., by invoking lowest-lying meson dominance and/or quark counting rules. 
In this context, assuming factorization, different phenomenological g^-dependences have 
been scrutinized in Ref. [^, attempting to reproduce the available data. 


Clearly, experimental studies of other color-suppressed B decays should be very desirable 

in order to complete the present information and more extensively test the factorization 

^For a rather complete set of references on heavy-to-light meson semileptonic form factors see, e.g., 
Refs. US, |. 

^Arguments in support of factorization for finite Nc, in the limit of infinitely heavy quarks, can be 
found in (with some criticisms in j^), and different analyses of the role of nonfactorizable matrix 
elements for selected B and D nonleptonic decays [|^, do not seem to indicate sizable deviations if 
1 /Nc corrections are collectively discarded [|^ . 
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scheme. From this point of view, an interesting role is played by the transitions B —> 
K(*) rjc, that differ from the previous ones, B K* J/'h, just by the charm-quark 
spin conhguration. These processes, not measured yet, have been considered also in the 
phenomenological analyses of Refs. 10. H. where relations between the branching ratios of 
B —> K^*'^ rjc and B K^*'> J/T in terms of ratios of leptonic constants and form factors 
have been worked out. 


In this note, we consider predictions for B —> K^*'^ r]c and B —> rj'^ that can be 

obtained in the factorization approximation using the method of QCD sum rules. The aim 
of our discussion differs from Bi in that we adopt here a specihc theoretical framework to 
estimate consistently from the outset all nonperturbative quantities necessary to predict 
these processes. Indeed, in this case the relevant leptonic constants in Eq. (^, where 
J/T must be replaced by rjc or rj'^, are 

< Vc (hc)|c7M75c|0 >= (^/) , (6) 


that are not experimentally known, but must be determined theoretically. Previous ref¬ 
erences rely on constituent quark model arguments to relate fr,^ to fj/q,. In fact, leptonic 


constants represent a natural field of application for QCD sum rules [^. The virtue 
of this approach, being fully relativistic and held-theoretic by construction, is that it al¬ 
lows to avoid the notion of conhned-quark wavefunction, and it incorporates naturally 
fundamental features of QCD, such as perturbative asymptotic freedom and nonpertur¬ 
bative quark and gluon condensation. Thus, as a computational scheme alternative to 
the constituent quark model, it should be sensible to present an updated estimate of frf^ 
and /,,/ in the framework of QCD sum rules. We recall that besides their relevance to 
the B decays under consideration, emphasized in Ref. |]^|, these constants are of phe¬ 
nomenological interest also to other applications of QCD sum rules, such as the radiative 
processes J/T —> ?7c 7 and rjc {j}'^ —>■ 77 H- 


The other signihcant aspect of the present calculation is that the needed B —> K^*'^ form 
factors Foi^q^) and Aoi^q^) (see Eqs. (^) and (H)) have not been studied as extensively as 
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others. Since these form factors must be evaluated at or we can exploit 

the capability of QCD sum rules to provide an evaluation of their g^-dependence, within 
some theoretical uncertainty. This would complete the phenomenological description of 
color-suppressed B decays to charmonium based on QCD sum rules and factorization, 
by the addition of a set of predictions which, ultimately, could be interesting to compare 
to experimental data (when available) or to predictions from alternative nonperturbative 
schemes. 

The paper is organized as follows. In Sec. 2 we present a QCD sum rule calculation of the 
leptonic constants and together with a comparison with previous determinations. 
In Sec. 3 we briefly describe a calculation of the form factors To(g^) and for 

the transitions B ^ K, K*] finally, in Sec. 4 we give predictions for ratios such as 
T{B K*r]'^)/T{B K*t]^) and r(5 ^ Kr]'J/T{B Kr]^). 


2. QCD sum rule calculation of L and /„ 


In order to determine the leptonic constants of the rjc and t]), mesons from QCD sum rules 
we consider the two-point function 

Aiq) =ijd^x < 0|T(aM^(a;)aMt(0))|0 > , (7) 

where d^A^{x) = 2mc : c{x)i'^^c{x) :. This two-point function is known in perturbative 
QCD to two-loop order |^, and the leading non-perturbative term in the Operator 
Product Expansion (OPE) is also known |]^ (it is given in terms of the gluon condensate). 
We exploit two different types of QCD sum rules, viz. Hilbert transforms at = 0, and 
Laplace transforms. The former can be written as 

1 ds 


rOO rja 

= ;7T7l^V'5(s) . 


The perturbative QCD contribution is given by 

3 r2(n) 1 


71—1 


M^(0)|p^ — 


r(2n) (2n -|- 1) \m‘l 


- - (1 -f- ttnOis) , 


( 8 ) 


(9) 
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where 




Svr 


1 - 


1 1 
+ ( - + 


2n n + 1 2(n + 2) \n n + 1 


n+1 

xE 

r=l ^ 


n + 2 B{r,\) _ 3 ^ ! 

r B(n, \) r 2r 


( 10 ) 


and B [x, y) is the beta function. The leading non-perturbative term in the OPE gives 

(n-3)r(i)r(n + 3) 


the following contribution |16 


MniO)\NP 24 ^ TT ^ (4m2)-+i r(n + 5/2) 
The hadronic spectral function is parametrized by the standard ansatz 

^ Im ^l^5{s)\HAD = (s - MjJ + e{s - So) ^ Im ^(’5(5)I 


( 11 ) 


TT 


PT 


( 12 ) 


where a second pole term for the t]'^ may be added to Eq. ([T^). Putting everything to¬ 
gether, and ignoring as a hrst step the contribution of the 77 ' meson, the Hilbert transform 
QCD sum rules become 

f^o ds 1 ^ 1 as ^2 1 


f2 
> ric 




X 


(n-3)P(l)P(n + 3) 

r(n + 5/2) 


(13) 


We take ratios of two consecutive moments for various values of n in order to hnd the 
duality window (range of values of Sq) for which the mass of the rjc is correctly predicted. 
Once this mass is accounted for, the leptonic decay constant follows from any given 
moment. In addition to expecting some region of stability in sq, one also expects the 
mass and decay constant to be reasonably independent of n. These expectations are 
partly fulhlled by the results of our calculation, as may be appreciated from Figs. 1 and 
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2 where we show Mrj^ and = frj^/V^ as a function of n for the particular choice: 
A = 200 MeV, rric = 1.39 GeV, and < >= 0.063 GeV"^. The solid (dash) curve in 

Fig. 1 is the theoretical (experimental) result. In Fig. 2 the solid (dash) curve is obtained 
by using the theoretical (experimental) value of the t]c mass in the Hilbert moments. 
While the stability in n is more than adequate, there is some sensitivity to changes in 
sq. This feature is well known from applications of QGD sum rules to heavy-light quark 
systems. The situation here, in the case of two heavy quarks, is then not much different. 
We have explored the following range of parameters, dictated by the gluon condensate 
and quark-mass analyses of pT|, respectively: < >= 0.063 — 0.19 GeV^, 

rric = 1.46 ± 0.07 GeV. The QGD scale was varied in the range A = 200 — 300 MeV. 
The values of sq needed to correctly reproduce the experimental value of were found 
in the range Sq = 9.5 — 10.5 GeV^. Results for and corresponding to various 
choices of A, rric, the gluon condensate and sq are qualitatively similar to Figs. 1 and 2. 
Quantitatively, we hnd frj^ ~ 256 — 300 MeV. 


In a second step, we include a second pole in the spectral function Eq. (|^ to account 
for the 77 ' meson, and repeat the analysis. We find that, in order to reproduce correctly 
the experimental values of both and (for the mass of r;' we use the value = 
3595 ± 5 MeV [^]), the value of Sq must increase by some 30% (sq ~ 13 — 14 GeV^). 
This increase is to be expected, since now Sq > Mi. We find: 

'Ic 

fr,^ ~ 301 - 326 MeV, ~ 231 - 255 MeV (14) 

and, therefore, the result for is now somewhat higher. The reason why the uncertainty 
in is now smaller than before is easy to understand. With a two-pole spectral function, 
the sum rules are now constrained to reproduce correctly the two meson masses simulta¬ 
neously. This constraint effectively reduces the size of the parameter space spanned by 
A, rUc, < > and Sq. 
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Next, we make use of the Laplace transform QCD sum rule and its first derivative ^3 

= f ds - ImV's(s) IPT 

J4m'i TT 


6y/7l 




3 3 4mr 


2 ’ 2 ’ Ml 


-6G 


1 3 4m?' 


2 ’ 2 ’ Ml 


(15) 


= r dsse-^/^l-lm'il,{s)\pT 

J Arrii TT 


A 2 ^2 -4m2/M2 3 3 4m?' 

-4m: < ctsG^ > e G --, -, 


6i/? 


2 ’ 2 ’ Ml 


.6G (-11 Li'l - Ig (-11 Li'l + 3G (1. ?. 


2’ 2’ M|y 2 I 2’ 2’ Ml 


2' 2’ Ml 


(16) 


where Ml is the Laplace variable and the function G{h^ c, u) 


G{b, c, uj) = 


00 


/ 1 + - 
i (cj Jo \ 00 


-b 


(17) 


is related to the Whittaker function 1 La^^(ci;) through 

G{h,c,u) = w'-i . 


(18) 


Considering hrst a one-pole spectral function, the ratio of Eqs. (|T^) and (pISD determines 
Mrj^ as a function of the Laplace variable Ml and of Sq. We hnd that is stable against 
changes in Ml in the wide region Ml ~ 1 — 5 GeV^. The experimental value of M^^^ is 
correctly reproduced for sq — 10 GeV^, roughly independent of the values of A, rric and 
< asG^ >. This value of sq is in agreement with the one obtained from the Hilbert 
moments. In Fig. 3 we show the result for (solid line) compared to the experimental 
value (dash line), for A = 200 MeV, rUc = 1.39 GeV and < agG^ >= 0.063 GeV"^. The 
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result for obtained by using the predicted (experimental) value of is 

shown in Fig. 4 as the solid (dash) curve. Exploring the parameter space of A, rric, 
< > and sq leads to the prediction: ~ 265 — 274 MeV, in agreement with the 

Hilbert moments. 


Adding a second pole to the hadronic spectral function, considering higher derivatives of 
Eq. (|T^ in 1/Ml, redoing the analysis we hnd that it is possible to reproduce the 
experimental values of both and ikh,/ for sq — 13 — 14 GeV^. At the same time, 
the stability region in the Laplace variable remains wide: ~ 2 — 5 GeV^. The 

predictions for the leptonic decay constants are 

~ 292 - 310 MeV, ~ 247 - 269 MeV , (19) 

in agreement with the results from the Hilbert moments, Eq. (14). Gombining the 
predictions from the two methods we obtain our hnal result 


fn, = 309 ± 17 MeV, /^/ = 250 ± 19 MeV , 

^ = 0.8 ± 0.1 . 

J'Hc 


( 20 ) 

( 21 ) 


Estimates of fn^, and of the J/T leptonic constant were made many years ago, in the 
early days of the QGD sum rules approach [^, |^. For example, in [^] a QGD sum 


rule determination of has been carried out at a one-loop approximation in perturbative 
QGD, and no gluon condensate corrections, in the framework of Hilbert moments. With 
an input value of the charm quark mass: rric = 1.25 GeV (to be compared with the recent 
determination employed here: rric = 1.46 ±0.07 GeV ), and the choice Sq = 16 GeV^, 
the value = 252 — 369 MeV has been obtained. Our determination here represents 
an improvement in many respects: radiative (two-loop) as well as non-perturbative cor¬ 
rections have been incorporated into the theoretical side of the sum rules, Hilbert and 
Laplace sum rules have been employed, updated values of the input parameters have 
been adopted, and the analysis was constrained to reproduce correctly the masses of both 
pseudoscalar mesons: and 






3 . FQ{q‘^) and ^0(9^) from QCD sum rules 


The form factors and AqIq"^), related to the weak transitions B ^ K and B K*, 

can be computed by three-point function QCD sum rules, following the same strategy 
adopted in the calculation of the leptonic constants fn^ and The starting point is 
provided by the correlators [^, ^ : 

n^uip,p',q) ^fjdx dy edr'-^-P-v) < 0|r{j„'''(i)l/„(0)jJ(!/)}|0 > (22) 

and 

Fu{p,p\q) =i‘^ J dx dy < 0\T{jI^* {x)A^{0)jl{y)}\0 > ; (23) 

the quark currents js and j^’^* are given by: js = dij^b, = dq^q^s, j^* = dj^s, 
whereas the flavour changing weak currents V and A are: = sq^b and A^ = s'y^q^b. 


The products: and receive a hadronic contribution from the states 5, K and 

B,K*, respectively, so that they can be expressed in terms of the form factors Fo{q‘^) 
and and of a continuum of states; on the other hand, a QCD calculation can be 

performed for the same operator products, in the limit of large and spacelike p^ and 
in terms of the perturbative QCD contribution and non-perturbative power corrections, 
proportional to vacuum matrix elements of high dimensional operators. The matching of 
the hadronic and QCD representations of the correlators can be improved by a double 
Borel transform in the variables —p^ and —We omit here the details of the calculation, 
and simply present the result of such a procedure for the form factor To(g^). Including 
power corrections in the OPE up to dimension 5 condensates, one obtains the following 
equation: 


fKfB — {Ml-Ml)Fo{q^) 

mb 




.){ 


mb — ms)\mb < qq > 


— j^dsds'p{s,s',q^) exp, 
ds < qgsF'^Gl^^q > } exp^^^ 


s Ml 
- + ^ 


M'2 


ml 


M2 


+ 


Ml 


mz 


M'2 


(24) 
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where the spectral function p is given by: 

1 r. . . (An - 2A's) 


p{s,s',q‘^) = ms)AH - - - 2A+2A'+n —2s')+ms(2s—n)]| (25) 

(A = s — ml, A' = s' — ml, u = s + s' — and \ = — 4 s s'). The coefficient reads: 


dfi — 


mlmg {mb — mg) 2mb mf mb{ml + ml — mbmg — q^) 


+ 


+ 


+ 


(26) 


4M'4 6M'2 3M2 4M4 

The integration region D in Eq. (|^) is limited by the effective thresholds sq and Sq 
separating, in the sum rule, the contribution of the ground state from the hadronic con¬ 
tinuum, the latter being modeled by the leading order perturbative QCD. The param¬ 
eters employed in the calculation are the quark masses: mb = 4.6 GeV [0 and = 
0.175 ± 0.020 GeV (we neglect up and down quark masses), the leptonic constants 
Jb = 0.18 GeV 1^ and fx = 0.16 GeV. As for the dimension 3 and dimension 5 con¬ 
densates, we take < gg > (1 GeV) = (—230 MeV)^ and < qgsO''"'G1^,^^q >= m^ < qq >, 
with tUq = 0.8 GeV^] within the hnal uncertainties, rescaling the quark and mixed con¬ 
densate to higher scales by the leading-log approximation of the anomalous dimensions 
does not affect the numerical results for the form factors. 


To derive Fq from Eq. (plf) , one looks for a region where the result is stable under variation 
of the Borel parameters M^, M"^, and of the continuum thresholds Sq for the B-cha nn el 
and Sg for the TC-channel. We hnd that stability is obtained by choosing the continuum 
thresholds in the ranges Sq = 33 — 36 GeV"^, Sg = 1.3 —1.5 GeV^, and the Borel parameters 
in the range = 8 ± 1 GeV^ and M"^ = 2.0 ± 0.4 GeV^. Let us hnally notice that the 
OPE expansion of the correlators in Eqs. (P^D and (P^D, starting from spacelike values 
of the momentum transferred, can be extrapolated to positive g^ provided that one is far 
from non-Landau singularities We compute the form factors up to g^ = 15 GeV^, 
which is in the safe region. 

The result for Fo{q‘^) corresponding to the above input numbers is depicted in Fig. 5. 


The analogous expressions relevant to Ag(g^) can be found in Ref. |Q. In this case, the 
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same input parameters are used as for Fo(g^), except that we need, in this case, fx* = 
0.22 GeV and the threshold in the i^*-channel is now in the range Sq = 1.5 — 1.7 GeV^ 
(while the range of Sq for the f?-cha nn el is the same as before). The results for Ao(q'^) are 
displayed in Fig. 6. 


As an important virtue of this calculation, we remark that for both Aq and Fq there 
exists, in the sum rules, a hierarchical structure in the OPE expansion, namely the leading 
perturbative term is numerically larger than the power corrections. This is a particularly 
welcomed feature of the present calculation, from the point of view of this theoretical 
framework, in comparison e.g. with the analogous determination of the form factor A 2 {q^). 


Turning to a discussion of the results, we observe in AQ{q^) a sharp increase when varies 
in the range = 0 — 15 GeV^. The form factor can be fitted by the expression: 


Mq^) = 


^( 0 ) 

X _aL 

^ Ml 


(27) 


with Ao(0) = 0.27 ± 0.03 and Mp = 4.8 ± 0.2 GeV. Therefore, the mass of the pole 
in Aq is slightly smaller than the mass of the hrst singularity in the q^ channel: = 


5.375 ± 0.006 GeV. This behaviour is similar to what has been observed in Ref. [Q for 
the channels D —> K*, D —y p, and B ^ p. 


On the other hand, the g^-dependence of F^^q^) is rather soft. We obtain Fo(0) = 0.29 ± 
0.03, whereas the fitted mass of the pole in Eq. (0) is Mp ~ 7.5 GeV, to be compared to 
the mass of the first resonance in the g^ channel, a bs{0^) state expected, in constituent 
quark models, in the region near 6 GeV (in the BWS model the value: M(fes)(0^) = 
5.89 GeV is chosen for the mass of this state). This behaviour is analogous to that 
observed for the form factor Fo{q‘^) in B —y 7i, computed in the framework of the inhnite 
heavy quark mass limit in Ref. . 


These results support the observation, made in Ref. |0, that the assumption of a given 
g^-dependence (such as polar, multipolar, etc.) in heavy meson semileptonic form factors 
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should be more safely confirmed by some explicit theoretical calculation. In fact, although 
the dominance of a singularity is reasonable for high values of q^, it is possible that, for 
lower q^, the contribution of higher states or of different dynamical mechanisms modifies 
the polar behaviour. The same argument holds for the application of QCD counting rules, 
that a priori are rigorous in the limit of large, spacelike 


4. Predictions and conclusions 


In the constituent quark model the leptonic constants of the charmonium system can be 
expressed in terms of the cc wave function at the origin T(0) ||17|| : 


TY)^ 

/l=48^l^(0)r 

(28) 



Mj/^ 

(29) 


therefore, the ratio fr^^/fj/m can be predicted in terms of the meson masses and of the 
charm quark mass: 


^ = =0.97 ±0.03, (30) 

fj/^ V Ml ^ 

where the value for rric chosen in Sec. (2) has been used. It is known that relativistic and 
radiative QCD corrections could modify the prediction in Eq. (|30D. Using our result in 
Eq. (PPD, which includes both, and the experimental value fj/q, = 384 ± 14 MeV, we get: 


f 

^ = 0.81 ±0.05. (31) 

Jj/i) 

A comparison of this result with that in Eq. (|30[) suggests that corrections are at the 
level of 15 — 20%. On the other hand, Eq. (^), applied to the radial excitations rj'^ and 
4/', gives the prediction: ^ = 0.80 ± 0.02, to be compared to ^ = 0.87 ± 0.08, obtained 
from Eq. (p0|) and using the experimental measurement f^' = 282 ± 14 MeV, possibly 
suggesting a minor role of radiative and relativistic corrections for these states. 


Turning to nonleptonic B decays, we first analyze ratios of decay widths, such as B ^ 

KV) 

T]c and B in the factorization approximation, where the dependence on 
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the Wilson coefficient 02 and on other weak parameters drops out, the relevant remaining 
dynamical quantities being the leptonic constants and the semileptonic form factors. 


Consider, for example, the ratio: 

r(5- - 


Rk = 


K-,i) 


= 0.771 f^l 




(32) 


l?c7 


where the numerical term is a phase space factor. The interesting point is that, because 
of the flat shape of id)(g^), Rk mainly depends on the ratio of the leptonic constants: 


Rk = 0.771 (1.09 ± 0.09)^ = 0.60 ± 0.15 . 

Jr)c 


(33) 


Thus, in the factorization approximation, a measurement of Rk would provide us with 
interesting information on and complement our knowledge of the properties of the qq 

J'Hc ‘ 

wavefunction. 


On the other hand, the analogous ratio for the decays into K* is given by 

I I n —> t\ rr \ r„/ r. j^n\ivR 

Rk* = 


r(B- ^ K-4) _ „ j,... 


r{B-^K-r,,) b„(Mp 


)= 0.381 (1.4 ±0.2)W (34) 

Jrjc 


Here, the ratio of the form factors deviates from unity due to the g^-dependence of Aq. 
The prediction from (^) would be: Rk* = 0.45 ± 0.16. Moreover, we observe that 


in the assumed factorization approximation the quantity \JRk*/Rk is sensitive to the 
g^-dependence of the ratio Aq/Fq. 


1.42 


Rk* 

,7l„(M.|)/fi,(M4). 

\ Rk 



(35) 


i.e. mainly to the g^-dependence of Aq since Fq is predicted to have a rather flat behaviour. 


A bound on the ratio 


R 

r(H- 




= 0.373 


K-7],) ^Fo(M,2 


(36) 


l?c7 


has been derived in Ref. |@|: 0.19 < Rr^^ < 0.98. This bound is satisfied by our result R^/^ = 
0.73 ± 0.13. In addition, for the analogous quantity 7?,^/, we predict: = 0.56 ± 0.12. 
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Finally, we consider the ratio: 


r(5- ^ k-7^,) 
T{B- K-J/^) 


2.519 

Jj/'S> 


mil) '2 


(37) 


In order to predict Rx we can use the relation: Fi(0) = Fo(0), and the observation, 
common to the QCD sum rules calculations of the g^-dependence of Fi, of the validity 
of the single pole model dominated by the 1“ B* resonance [^, ^ 


361. We obtain: 


Rk = 0.94 ± 0.25, and, for the analogous quantity R'x = • ^'k — 1-61 ± 0.53. 

This implies that, using the CLEOII experimental measurements: B{B~ K~ J/'^) = 
(0.11 ± 0.01 ± 0.01) X 10-2 and B{B- R-^') = (0.06 ± 0.02 ± 0.01) x lO’^ we expect: 
B{B- R-r]^) = (0.11 ± 0.03) x lO-^ and B{B- R-j]'^) = (0.10 ± 0.05) x lO'^, and 
therefore these decays are in a range well-accessible to the present experimental facilities. 


With these predictions we conclude our analysis. The relatively large decay rates of the 
processes B R^*^t]c and B R^*'>t]'^ should allow their observation in the near future. 
This measurement will shed more light on the problem of factorization, which is a basic 
assumption in the present analysis of heavy meson nonleptonic decays. 


From the theoretical point of view, it should be interesting to compare the results obtained 
here with other QCD calculations. For example, it should be possible to calculate the 
leptonic constants of t]c and t]'^ mesons by lattice QCD. Moreover, regarding the B —> R^*'^ 
form factors, a possible independent test of the g^.^ependence could be obtained, e.g., 
by using light-cone sum rules. Also, extrapolations of form factors in g^ and in the 


heavy quark mass, starting from the charm mass, in the lattice QCD framework |37 


could possibly provide us with enough information, at least at a qualitative level, to be 
compared with the QCD sum rules results. 


14 














Acknowledgments 

We would like to thank F.Buccella, F.De Fazio, G.Nardulli and P.Santorelli for interest¬ 
ing discussions. The work of (NP) was snpported in part by the HCM, EEC Contract 
ERBCHRXCT930132. The work of (CAD) was supported in part by the Foundation 
for Research Development (ZA) and the John Simon Cnggenheim Memorial Fonndation 
(USA). 


15 



References 


[1] M. Gourdin, A.N. Kamal and X.Y. Pham, Phys. Rev. Lett. 73 (1994) 3355. 

[2] R. Aleksan, A. Le Yaouanc, L. Oliver, O. Pene and J.C. Raynal, Report 

DAPNIA/SPP/94-24, LPTHE-Orsay 94/15. 

[3] For a review see: T.E. Browder, K. Honsheid and S. Playfer, in; ”B decays”, S. Stone 
Ed, World Scientific, 1994. 

[4] G. Altarelli, G. Gurci, G. Martinelli and S. Petrarca, Nncl. Phys. B 187 (1981) 461; 
A.J. Buras, M. Jamin, M.E. Lantenbacher and P.E. Weisz, Nncl.Phys. B 370 (1992) 
69 and Nncl.Phys. B 375 (1992) 501 (addendnm); 

M. Giuchini, E. Franco, G. Martinelli and L. Reina, Nncl. Phys. B 415 (1994) 403. 

[5] A.J. Buras, Report MPI-PhT/94-60, TUM-T31-75/94. 

[6] M. Wirbel, B. Stech and M. Bauer, Z. Phys. G 29 (1985) 637. 

[7] M. Gourdin, Y.Y. Keum and X.Y. Pham, Report PAR/LPTHE/95-01. 

[8] M. Gourdin, Y.Y. Keum and X.Y. Pham, Report PAR/LPTHE/94-32. 

[9] A.N. Kamal and A.B. Santra, Report Alberta Thy-27-94 and Alberta Thy-31-94; 

G. E. Garlson and J. Milana, Report WM-94-110; 

H. -Y. Gheng and B.Tseng, Report IP-ASTP-21-94. 

[10] G. ’t Hooft, Nncl. Phys. B 72 (1974) 461; 

E. Witten, Nncl. Phys. B 149 (1979) 285. 

[11] M. Dugan and B. Grinstein, Phys. Lett. B 255 (1991) 583. 

[12] U. Aglietti, Phys. Lett. B 292 (1992) 424. 

[13] B. Block and M.A. Shifman, Sov. Jour. Nncl. Phys. 45 (1987) 135, 301, 522; Nncl. 
Phys. B 389 (1993) 534, Nncl. Phys. B399 (1993) 441, 459. 


16 



[14] A. Khodjamirian and R. Riickl, Report MPI-PhT/94-26, LMU 05/94. 

[15] M. Bauer, B. Stech and M. Wirbel, Z. Phys. C 34 (1987) 103. 

[16] M.A. Shifman, A.I. Vainshtein and V.I. Zakharov, Nucl. Phys. B 147 (1979) 385,448; 
L.J. Reinders, H. Rubinstein and S. Yazaki, Phys. Rep. 127 (1985) 1. 

[17] N.G. Deshpande and J. Trampetic, Phys. Lett. B 339 (1994) 270. 

[18] V.A. Beilin and A.V. Radyushkin, Nucl. Phys. B 260 (1985) 61. 

[19] D. J. Broadhurst and S. C. Generalis, Open University Report No. OUT-4201-8 
(1982); 

S. G. Generalis, PhD thesis OUT-4102-13 (Open University, 1983); 

see also: K. Schilcher, M.D. Tran and N. F. Nasrallah, Nucl. Phys. B 181 (1981) 91; 

B 187 (1981) 594 (E). 

[20] R.A. Bertlmann et al. Z. Phys. G 39 (1988) 231; 

G.A. Dominguez and J. Sola, Z. Phys. G 40 (1988) 63. 

[21] G.A. Dominguez, G.R. Gluckman and N. Paver, Phys. Lett. B 333 (1994) 184; 

S. Narison, Report GERN/TH-7444/94, to appear in the Proceedings of the 
”QGD94” Workshop, Montpellier (France) 7-13 July 1994. 

[22] G. Edwards et al., Phys. Rev. Lett. 48 (1982) 70. 

[23] R.A. Bertlmann, Nucl. Phys. B 204 (1982) 387. 

[24] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions (Dover, New 
York, 1972). 

[25] V.A. Novikov et al, Phys. Rev. Lett. 38 (1977) 626; Phys. Lett. B 67 (1977) 409; 
Phys. Rep. 41 (1978) 1. 

[26] L.J. Reinders, H. Rubinstein and S. Yazaki, Phys. Lett. B 113 (1982) 411. 


17 



[27] M. Shifman, JETP Lett. 30 (1979) 546. 

[28] B.L. Ioffe and A.V. Smilga, Nucl. Phys. B 216 (1983) 373; 

A.V. Nesterenko and A.V. Radyushkin, Phys. Lett. B 115 (1982) 410. 

[29] For a review see the reprints volume: ’’Vacuum structure and QCD Sum Rules”, 
edited by M.A. Shifman, North-Holland, 1992. 

[30] C.A.Dominguez and N.Paver, Phys. Lett. B 293 (1992) 197. 

[31] C.A.Dominguez, C. van Gend and N.Paver, Phys. Lett. B 253 (1991) 241; 
K.G.Chetyrkin, G.A.Dominguez, D.Pirjol and K.Schilcher, Report MZ-TH/94-21 
(1994); 

M.Jamin and M.Miinz, Report GERN-TH.9435/94 (1994). 

[32] G.A.Dominguez, in Proceedings of the Third Workshop on the Tau-Gharm Factory, 
Marbella, Spain, 1993, edited by J.Kirkby and R.Kirkby, Editions Frontieres, Gif- 
sur-Yvette, France, 1994. 

[33] P. Ball, V. Braun and H.G. Dosch, Phys. Rev. D 44 (1991) 3567. 

[34] P. Golangelo, F. De Fazio and P. Santorelli, Report BARI-TH/94-174, to appear in 
Phys.Rev.D. 

[35] P. Golangelo and P. Santorelli, Phys. Lett. B 327 (1994) 123. 

[36] P. Ball, Phys. Rev. D 48 (1993) 3190; 

V.M. Belyaev, A. Kodjamirian and R. Riickl, Z.Phys G 60 (1993) 349. 

[37] A.Abada et ah, Nucl. Phys. B416 (1994) 675; 

K.G.Bowler et ah, UKQGD Gollaboration, Edinburgh Preprint 94/546, Southampton 
Preprint SHEP 93/94-32. 


18 



FIGURE CAPTIONS 


Fig. 1 

The mass of the rjc from Hilbert moments as a function of n for A = 200 MeV, rric = 
1.39 GeV, and < >= 0.063 GeV"^. The solid curve is the prediction and the dash 

curve the experimental value. 


Fig. 2 

The leptonic decay constant using a single pole spectral function in Hilbert 

moments, as a function of n and for the same values of input parameters as in Fig. 1. 
The solid (dash) curve is obtained using the predicted (experimental) value of 


Fig. 3 

The mass of the rjc from the Laplace transform for A = 200 MeV, rric = 1.39 GeV, 
and < >= 0.063 GeV^. The solid curve is the prediction and the dash curve the 

experimental value. 


Fig. 4 

The leptonic decay constant fn^ = fr^^/^/2, using a single pole spectral function in the 
Laplace transform, and for the same values of input parameters as in Fig. 3. The solid 
(dash) curve is obtained using the predicted (experimental) value of Mr/^. 
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Fig. 5 


The form factor -Fo(g^) for the transition B ^ K. The curves refer to the sets of 
parameters: Sq = 33 GeV^ and Sg = 1.3 GeV'^ (continuous line), Sq = 33 GeV^ 
and Sg = 1.5 GeV^ (dashed line), sq = 36 GeV^ and Sg = 1.3 GeV^ (dotted line). 
So = 36 GeV^ and Sg = 1.5 GeV^ (dashed-dotted line). The Borel parameters in the B 
and K channel are fixed to = 8 GeV^, = 2 GeV^, respectively. 


Fig. 6 

The form factor AQ{q^) for the transition B K*. The curves refer to the sets of 
parameters: sq = 33 GeV^ and Sg = 1.5 GeV^ (continuous line), sq = 33 GeV^ and Sg = 
1.7 GeV^ (dashed line), sq = 36 GeV^ and Sg = 1.5 GeV^ (dotted line), sq = 36 GeV"^ 
and Sg = 1.7 GeV^ (dashed-dotted line). The Borel parameters in the B and K* channel 
are the same as in Fig. 5. 


20 



This figure "fig 1-1.png" is available in "png" format from: 


http://arXiv.org/ps/hep-ph/9502374vl 



This figure "fig 1-2.png" is available in "png" format from: 


http://arXiv.org/ps/hep-ph/9502374vl 



This figure "fig 1-3.png" is available in "png" format from: 


http://arXiv.org/ps/hep-ph/9502374vl 



